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Abstract 


We consider a fully-discrete approximation of the Allen—Cahn equation, 
such that the forward Euler/Crank—Nicolson scheme (in time) combined with 
the RBF collocation method based on “shifted” surface spline (in space). 
Numerical solvability and stability of the method, by using second order finite 
difference matrices are discussed. We show that, in the proposed scheme, the 
nonlinear term can be treated explicitly and the resultant numerical scheme is 
linear and easy to implement. Numerical results that show the efficiency and 
reliability of the proposed method are presented, and two types of collocation 
nodes for solving this equation are compared. 
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1 Introduction 


The radial basis function method is one of the most useful tools for interpolat- 
ing multidimensional scattered data. The method has remarkable properties, 
to handle arbitrarily scattered data, to adjust to several space dimensions 
easily, and to provide powerful convergence properties, which have made it 
popular in sciences and mathematics; see [16, 26, 6, 28]. Some of the more 
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recent applications of radial basis functions include measurements of poten- 
tial or temperature on the earth’s surface [3], mathematical finance [11], data 
mining [19], neural network [18], tomographic reconstruction, medical imag- 
ing [4], and the numerical solution of partial differential equations (PDEs). 


The Allen—Cahn equation was initially introduced by Allen and Cahn [1] 
in 1979 for antiphase domain coarsening in a binary alloy. This equation, as 
the simplest model of phase-field, describes the evaluation of diffuse phase 
boundary, concentrated in a small region. It is known that the Allen—Cahn 
equation is well-posed and has a unique solution that satisfies the maximum 
principle; see [12]. The Allen—Cahn equation has been widely used in many 
complicated moving interface problems in fluid dynamics and material science 
applications, such as crystal growth [7] and tumor growth [27] on flat surfaces. 


During the past two decades, a great deal of mathematical efforts have 
been devoted to the study of numerical solution for the Allen—Cahn equation 
[20, 31, 12], and finite difference [13], finite element [9], Fourier spectral 
with periodic boundary condition [15], and RBF methods[17] are some of the 
recently proposed methods. 


In this paper, a fully-discrete approximation of the Allen—Cahn equation 
is considered, such that a first order (Euler/Crank—Nicolson) finite difference 
scheme (in time) combined with the RBF collocation method (in space). 
Here, the main goal is how to construct an efficient numerical scheme based 
on the RBF collocation method that preserves the maximum principle and 
energy dissipation law in the Allen—Cahn equation. We will use a special kind 
of RBFs as the shifted surface splines [8, 29, 28, 30]. These functions include 
all of the multiquadric, thin plate, and power RBFs, which their properties 
quite well understood, theoretically and practically [3, 8, 26]. For example, 
some features of multiquadric, thin plate, and power RBFs are as follows: 

e Multiquadrics are strictly positive definite, but surface splines and powers 
are conditionally positive definite. 

e Powers and surface splines give an algebraic rate of convergence, but mul- 
tiquadrics can achieve spectral rate of convergence. 

e Multiquadrics interpolation is always solvable, provided that the points are 
disjoint, but for surface splines and powers, for solvability, we need to add 
some appropriate number of functions (polynomials). It is worth noting the 
shifted surface spline (—1)"~*/?++ (a? + gyre? log (x? + 6?) \/? where s 
is even and 6 ¥ 0, is strictly conditionally positive definite (see [8, p. 104]). 
Hence, for the well-posedness of the interpolation problem, we need to add 
some appropriate polynomial to this function. 

e Powers and surface splines are piecewise smooth, but multiqudrics are in- 
finitely smooth. 

e Franke’s research [10] showed that multiquadrics, in terms of simplic- 
ity of implementation, stability, and accuracy, have much better perfor- 
mance among 29 available interpolation methods used for interpolating two- 
dimensional analytic functions. 
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In the proposed scheme, the nonlinear term can be treated explicitly such 
that, the resultant numerical scheme will be linear and easy to implement. 
The numerical solvability and stability of the method by using the RBF-FD 
procedure based on the classical finite difference method are discussed. The 
RBF-FD method uses the simplicity of FDM formulation and the meshfree 
property of RBFs, which was proposed in [22, 21]. Also, Tolstykh, Lipavskii, 
and Shirobokov [24] and Wright and Fornberg [25] introduced this method 
independently for the numerical approximation of the ODEs and PDEs. For 
flat RBF, Wright and Fornberg [25] showed that the RBF-FD scheme is the 
same as the conventional finite difference scheme. The RBF-FD method is 
used with other names such as local radial basis function method [23] and 
local radial basis functions based on the differential quadrature method. By 
using some numerical experiments, we show that the proposed method with 
equally spaced collocation nodes may give better results in comparison with 
scattered nodes. 


The paper is organized as follows. In Section 2, the RBF interpolation 
is introduced in brief. Section 3 reviews the RBF collocation methods. In 
Section 4, numerical results that show the efficiency of the proposed method 
are given. The last section ends with concluding remarks. 


2 The RBF interpolation method 
Before we pose the RBF interpolation, we need to recall the following defini- 
tions. 


Definition 1. A continuous and even function f : R? > R is said to be 
strictly conditionally positive definite of order k (SCPD(k)) if, for all N EN, 
all sets y = {Xo,..., Xn} C R4, and all A € R™ satisfying 


N 
i=0 
for all real-valued polynomials q of degree less than k, the quadratic form 
N N 
MAA= STS 0 AAS (Xi — X5ll), 
i=0 j=0 


is positive, unless is zero. 


Let Il, (R2) be the space of d-variate polynomials of degree at most k; 
then we have the following definition. 


Definition 2. A set of distinct points {Xo,...,X} in R® is said to be 
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Il, (IR¢)-unisolvent if the only element of I; (R*) that vanishes at each x; is 
the zero polynomial. 
The definition of multivariate RBF interpolation can be stated as follows: 


Definition 3. Given a set x = {Xi}, of distinct centers in R? and a 
target function v : R? > R, the RBF interpolant is given by 


N M 
Sv(X) = $0 Aso (IX — Xl) + Do rere (X), MSN, (1) 
j=0 k=0 


where ¢ (||X — X,]|) is a strictly conditionally positive definite RBF of order 
m, centered on the points z,;, and {p;, (X es forms a basis for the polyno- 


mial space II,,-1 (R4), where M = cert The expansion coefficient .,, 


‘Yr are determined by the interpolation conditions, 
Su(Xj) =v(Xj), OSISN, 


and the following additional conditions are considered to handle the addi- 
tional degree of freedom 


N 
So Apr (Xj) =0, O<Sk<M. 
j=0 


These conditions lead to a system of equations that can be written down in 
the matrix form as 
® P A\_[(V (2) 
Pt 0 y} \O}?’ 


Y= (90,---,Y¥u) , V=(vo,---, UN) - 


where 


As it is shown in [26] that system (2) is uniquely solvable when ¢ is strictly 
conditionally positive definite of order m and y is I, -1 (IR*)-unisolvent set 
of centers; so the coefficient matrix in (2) is invertible. 

In this paper, we use “shifted” surface splines as radial basis functions 
@, which include the multiquadric and thin plate spline and are defined as 
follows: 


ae (—1)ln-s/21 (a? + gyms? sis odd, 
i) = > m—s 
eae + 6?) f2 log (a? + 52)? sis even, 


where s,m € N, m > s/2 and [s] indicates the smallest integer greater than 
s. For odd s, ¢5 is called the multiquadric RBF, and for 6 = 0, the function 
@s5 is called the surface spline. 
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In what follows, we consider the number hf as the “density” of y = 
{Xo,...,Xw} in Q, as follows: 


h=h(Q;x) := sup min |X — X;|. 
XEQ XKEXx 


3 The RBF collocation method, stability and numerical 
solvability 


In this section, we use the forward Euler/Crank—Nicolson scheme for the time 
and classical RBF method for spatial discretization in numerical approxima- 
tion of the Allen—Cahn equation. Let us define t, = k *7,k = 0,1,...,M, 
where T = 5 and M is the number of the time steps. we introduce notations 
v;(X) = v(X,t;) and jv; = “+42. Also we define 0; = Mitr Ms 


Consider the following time-dependent Allen—Cahn equation: 


n 


a) _ Ay(X,t)+45f(v(X,))=0, XEQCR4 te (0,7) is 
v (X,t) = 9(X), X € AQ, 


which is supplemented with the following initial condition: 
vu (X,0) =v (X) ind, 


where 2 is a bounded convex domain with Lipschitz boundary 00, T > 0 
is a positive constant, and f(v) = v? — v is the reaction term with F(v) = 
4(v? — 1)?, a given energy potential that drives the solution into the two 
pure states v = +1. The nonlinearity f(v) is monotone outside the interval 
[—1, 1], so solutions of the Allen—Cahn equation satisfy a maximum principle; 
see [12, 2]. Also, g(X) and X € OQ are known constants. The parameter 1, 
which usually stands for the inter-facial width, is a small constant compared 
to the characteristic length of the laboratory scale. The function v (X,t) 
is a distribution function of the concentration for one of the two metallic 
components of the alloy. 


In the RBF collocation method, the solution v (X,t) is approximated by 
a linear combination of RBFs as 


N M 
w (X,t) = DOA OH O(WX — Xl) + Do re Ore (X), MEN, (4) 
j=0 k=0 


where @ and px(X) are defined in (1). The collocation points {X;}4_, has 
been selected in such a way that {Xo,..., Xy,} C Qand {Xy,41,.-.,Xw} C 
0Q. In RBF collocation method, the requirement is that (3) is satisfied 
exactly by (4) at the following relations. 
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t) — Any (Xist) + FoF (Ow (Xis8)) = 0. k =0,1,...,Nr, 
un (Xx, t) = g(Xx), k=N,+1,...,N, 

and the initial conditions are 

un (Xx,0) = v0, &=0,1,...,N. 


According to the forward Euler/Crank—Nicolson scheme in time, we dis- 
cretize the Allen—Cahn equation as follows: 


UN (Xx, t)41) — UN (Xx, t;) Aun (Xx, t)41) + Aun (Xx, t;) 
T 2 
1 
+ cad (uw (Xn, t3)) = 0, k=0,1,...,N, 


where 7 is a mesh size in time and X, = (aj, y:),i,1 = 0,...,n, with N = 
(n +1)? —1. Equivalently, 


(ow )ett — F(Avw)i** (om) + 5(Aow) 


— al (@w)k) =0, k=01,....N 6) 
where 
. ; j 
(vw)h = on (Xeyt), (Aon yh = Aon (Xeyt) = (Cow )an + On)yy)- 
In addition, we impose the following constraints, 
N 
S> Aj (t) pi (Xj) =0, OSES M. (6) 
j=0 


From (5) and (6), we have a system with N + M +2 unknowns {Ae} and 
{ye} and N+M-+2 equations. But, before imposing boundary conditions, 
we will simplify systems (5) and (6), by using finite difference matrices. 
Now, If we set ®j, = $(||X; — Xx||) and A®;, = A¢g(||X; — X;z||), then 
we can rewrite collocation equation (5) and (6) in the matrix notation as 


(28) 5) 5 2) 


+ (3 (ow) + (owt) . 


(7) 


: j+1 
where AJ+! = (45 J: Now, by using (2), we may write (7) as 
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(mo)-a(P VT) ) =e") 


where V! = [(uy)}, (un)i,.--, (uw) yl’, =0,...,M and 


® P 
A=(> ar 


which is a dense linear system and when solved, we obtain an RBF based on 
the numerical approximation for (3). However, for a stable computation, we 
have to consider, stability issues for time-dependent problems and adapt the 
proposed method for nonlinear equations. 


To obtain stability criterion for system (8), we will use a recently proposed 
local method that causes more flexibility for handling nonlinearities. In this 
method, essentially a classical finite difference (FD) method or generalization 
of them to scattered node layouts are considered, which is known as RBF- 
FD methods. This kind of methods have been investigated independently by 
Shu et al. [22, 21], Tolstykh, Lipavskii, and Shirobokov [24], Cecil, Qian, and 
Osher [5], Wright and Fornberg [25]. Here, we only consider its most gen- 
eral form, in which, the FD method consists of approximating second order 
derivative of base functions at a given point, based on a linear combination 
of the value of radial basis functions at some surrounding nodes. In other 
words, we will use a classical 2-dimensional case as follows: 


s (Xl) le= = Souitent (Xi), O<j<n, 1<k<n, 


where w(?) is called the FD weights and ¢, (||X|]) = $(||X — Xs||). This 
formula can be combined to create FD formulas for partial derivatives in two 
and higher dimensions. In classical central finite difference approximation of 
order two, FD weights can be considered as: 


k 


a4 (2) 
X||) |o=a, (|X; 
ages (IAI) le=w, 2 wi (|X|) 
_ (I Xp4all) — 20 (Xa) + bs (XG-aI) 
h? , 
where fh, is the spatial distance of nodes. In two-dimensional case, we have 
d? a bs (X51 ll) — 26s (IX5,all) + os (X5—-1, ll) 
7m” (|X|) |w=a; + ap (XID ly=ui _ ae = 


bs (I|X5,i-1ll) = 26s (XG all) + bs (54-11) 
h? ; 
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In the matrix notation, we may write A® = pF, where the finite 
difference matrix F() is as follows: 


-410...0 1 00 

1 00 1-4 1 0 
0 0 0.20.1 Os. 

2. . . . . . 
ae 0 00...0 0 00.. 
1 —41 00 10 
0 10...0 1 -41.. 


From A® = 4 F?)® and AP = 4 F)P, collocation equation (8) can be 
written as 


(m0) ae (fo 1) Ce o)J'( )-8 


where R is the right-hand side of (8). For example, for n = 2,(N = (n+1)?), 


we have 
-41 0 000 0 0 
Pea oe 
00010%1-41 0 
ALTO 1 60. Oo o 
F2=|/ 01041-41041 0 
0000%10i1-41 
0 1-41 1 00 0 
CO. 2 Oo +47 Oo 
coo 8 0 1 Oo 14 


Now, (9) can be written equivalently as 
[7 Ch ee Bt ee ae j)N 10 
N+1~ 553 | — af ((on)s) i 3 (un); + (uw )ele=o- (10) 


It is well known that the meee of FO) (A?) are \; = —4 sin” (sy); 


or p(F)) < 8, that is, | 1— 45 |< 1. Hence for a stable computation, we 
need a <1. 


Note that the eigenvalues of F(?) after i imposiie boundary conditions are 
less than or equal to # 0 for allO <i< N. Thus, 
clearly system (10) is numerically solvable. 


oy 


According to considered order of mesh points {2;, y;}, 


{{Zo, yot, {Zo, yi}, Les {z0, Yn}; {z1, yo}, {x1,y1}; oe {r1, Yn}; ae {In, Yn}} 
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and given boundary conditions on 29, 2n, Yo, Yn, system (8) after imposing 
boundary conditions is completed. This is a system with N + 1 unknowns 
and equations. Therefore, after imposing given boundary conditions, (8) can 


be solved directly, by using stability condition 75 < 1. 


4 Numerical illustrations 


In this section, we present four numerical examples to illustrate the accuracy 
and efficiency of the proposed scheme. According to our observations, when 
equally spaced nodes are used as a set of collocation points, we obtain better 
results in comparison with scattered mesh points. In all of test problems, we 
have used multiquadratic radial basis functions with tT < 7? and tT < h?/2. 
In Example 1, the standard surface spline is considered too. For comparison, 
at time t = T, the infinity norm of the error 


jv — Sv|lgg = max, jv (Xi, 7) — $v(XT)], 


and the root mean square (RMS) norm of the error 


1 N 


wD (w(K) - $0 (Xi, 7)”, 


i=1 


llv — Sully = 


are used, where N is the number of collocation points. 


Example 1.Consider two-dimensional Allen—Cahn equation with the fol- 
lowing initial and boundary conditions: 


Av + (v3 —v) =f, (x,y) € (-1,1)’, t> 0, 


v(+1,y,t) = 0, v(a#,+1,t) =0,t>0, 
v (x,y, 0) = 0, (x,y) € [=1, Nee 


where 7 = 100, the exact solution is u(x,y,t) = tanh ((x? — 1) (y? — 1) t) 
and the function f is computed accordingly. 


In this example, we apply the proposed method by using two types of 
nodes. As shown in Table 1, the obtained errors for equally spaced nodes 
are slightly better than scattered nodes. We plot the absolute error of this 
example in Figure 1 and the contour plot of the error in Figure 2. The 
comparison of errors of Example 1 by using different time steps for t = 10007 
at N = 100 is presented in Table 2. 
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Table 1: The comparison of errors of Example 1, by considering “shifted” surface spline, 
at t = 10007 by using scattered and equally spaced nodes where time step is T = 10~*. 


Equally Spaced Scattered 

N [lu — uv [e—anls [= uv] [e— “nla 

49 4.82043E — 6 1.41187E —5 4.26723E — 6 1.30836E — 5 
64 4.07495E — 6 1.33603.E — 5 3.35978E — 6 1.34205F — 5 
81 1.31830E — 6 5.20152E — 6 7.72547E — 7 3.28669E — 6 
100 1.12669E — 6 5.04211E — 6 7.65370E — 7 2.91935E — 6 
121 7.02509E — 7 3.32090E — 6 7.56430E — 7 1.94812E — 6 
144 6.73415E — 7 3.52046E — 6 6.22881 E — 7 2.15101E — 6 


Table 2: The comparison of errors of Example 1, by using various time steps at N = 100 
and t = 10007. 


Equally Spaced Scattered 
Time Step (7 =10-) | fu—unl [lw unly || lu — unl, Iu—unlp _| 
F 1.12669E — 6 5.04211E —6 || 7.65370E — 7 2.91935E — 6 
1/2x7 2.78156E — 7 1.17651E —6 || 1.33866E — 7 6.79208E — 7 
1/4x7 6.91386E — 8 2.85469E — 7 || 3.36881E — 8 1.65689E — 7 
1/8x7 1.72370E — 8 7.04178E —8 || 8.64618E — 9 4.10731E — 8 
1/16 x + 4.30344E — 9 1.749466 —8 || 2.19114F —9 1.02357E — 8 


Absolute Error 


Figure 1: The plot of the absolute error in Example 1, at t = 10007 and N = 100 with 
equally spaced nodes. 
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Figure 2: The contour plot of the error in Example 1, at t = 10007 and N = 100 with 
equally spaced nodes. 


In what follows, we solve Example 1, by considering the standard surface 


spline ¢(r) = r?logr. Comparing the results obtained in Tables 1-4, it 


can be seen that the “shifted” surface spline provides better results than the 
standard surface spline. The absolute error and contour plot of the error, for 
this case, is plotted in Figures 1 and 4, respectively. 


Table 3: The comparison of errors of Example 1, by considering standard surface spline, 


at t = 10007 by using scattered and equally spaced nodes where time step is T = 10~¢. 


Equally Spaced Scattered 

N [lu — uv [e— "nls [= ule [u— “nls 

49 6.51165E — 6 1.84565E — 5 5.57870E — 6 1.78168E — 5 
64 5.59294E — 6 1.74854 — 5 4.96732E — 6 1.79523E — 5 
81 1.84662E — 6 6.64194E — 6 1.13404E — 6 4.34977E — 6 
100 1.58833E — 6 6.32184E — 6 8.10437E — 7 3.54860E — 6 
121 7.04900E — 7 3.59848 E — 6 7.72250E — 7 1.96949E — 6 
144 6.74292 E — 7 3.73803E — 6 6.20674E — 7 2.14926E — 6 


Table 4: The comparison of errors of Example 1, by considering standard surface spline, 


by using various time steps at N = 100 and t = 10007. 


Equally Spaced Scattered 
Time Step (7 =10-) | __ u—unllax [u—unly | een [u-unly 
7 1.58833.E — 6 6.32184E — 6 || 8.10437E — 7 3.54860E — 6 
1/2x7 3.91406E — 7 1.50928F —6 || 1.85967E — 7 8.74747E — 7 
1/4xr 9.71936E — 8 3.69865E — 7 || 5.01409F — 8 2.20779E — 7 
1/8x 7 2.42193E — 8 9.16348E — 8 || 1.30555E — 8 5.57276E — 7 
1/16 x7 6.04515E — 9 2.28113E — 8 || 3.33334E —9 1.40169F — 8 
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Absolute Error 


- mt 
ly 


Figure 3: The plot of the absolute error in Example 1, by considering the standard 
surface spline, at t = 10007 and N = 100 with equally spaced nodes. 
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Figure 4: The contour plot of the error in Example 1, by considering the standard surface 
spline at t = 10007 and N = 100 with equally spaced nodes. 


Example 2. Consider the Allen—Cahn equation with the following initial 
and boundary conditions 


Ov Riscs 1 


3 _ 11) 
Ot ae v) bE (eye ( ’ ) ,t>0, 
v(£1,y,t) =e" sin ((y* — 1)t), v(w,£1,t) =0, ¢ > 0, 


v (x,y, 0) =0, (x,y) € [-1, 1] ; 


where 7 = 100, the exact solution is v(a,y,t) = e~™* sin ((y? — 1) t) and 
the right-hand side function f can be computed accordingly. 

The comparison of errors for two types of nodes is presented in Table 
5, and the plot of the absolute error is given in Figure 5. We also plot the 
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contour of the error in Figure 6. The obtained errors for this example, by 
using different time steps, are also presented in Table 6. 


Table 5: The comparison of errors of Example 2 for t = 10007 by using scattered and 


equally spaced nodes where time step is T = 10~¢. 


Equally Spaced Scattered 

N [he — Ung [u—unllp Te — unl [lu —unTlp 

49 1.07720F — 5 3.00790E — 5 1.31170F —5 3.15498E — 5 
64 8.74748 E — 6 2.65163E — 5 1.14132F —5 3.47482E — 5 
81 1.67367E — 6 4.61698E — 6 1.55723E — 6 4.72437E — 6 
100 1.33128E — 6 4.66584E — 6 1.43747E — 6 5.21502E — 6 
121 6.75815E — 7 2.56911E — 6 7.37549E — 7 1.53253E — 6 
144 6.28635E — 7 2.73066E — 6 5.62016E — 7 1.69418E — 6 


Table 6: The comparison of errors of Example 2 by using various time steps at N = 100 


and t = 10007. 
Equally Spaced Scattered 
Time Step (7 = 107) [u—unlln [u—unlp [Ie — anes [u—unlp _| 
a 1.33128E — 6 4.66584E — 6 || 1.43747E —6 5.21502E — 6 
1/2xr 3.29225E — 7 1.12282E —6 || 3.27643E —7 1.21081E — 6 
/4xr 8.18399E — 8 2.75750E — 7 || 7.82939F — 8 2.92148E — 7 
1/8x7 2.04010E — 8 6.83517E — 8 || 1.91410F — 8 T7.17766E — 8 
1/16 x 7 5.09283E — 9 1.70169F — 8 || 4.73236E —9 1.77901E — 8 
Absolute Error 

ae. ae ee) | 

ee i] 

| 

sso 
0.0001 VE a 
\ mE 


0.0000 


-1.0 


7 0.0 


Figure 5: The plot of the error in Example 2, at t = 10007 and N = 100 for equally 
spaced nodes. 
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Figure 6: The contour plot of the error in Example 2, at t = 10007 and N = 100 with 
equally spaced nodes. 


Example 3. For this example, consider the following equation: 


Ov 


v(+l1,y,t) =0, v(a,4 
v(x, y, 0) = sin (ry) cos (wx/2), (x,y) € [-1,1]?, 


H1,t)=0,t>0, 


5(v3—v) =f, (x,y) € (-1,1)*, t>0, 


(11) 


where 7 = 100 and the exact solution is v (x, y,t) = sin (7y) cos (7x/2) (1 + t) 
and the function f can be computed accordingly. 


The comparison of errors for the equally spaced and the scattered nodes 
for this example is shown in Table 7. As we expected, better results are 
obtained for equally spaced nodes. The plot of absolute error, in the case of 
equally spaced nodes, is presented in Figure 7, and the contour plot of the 
error is given in Figure 8. In Table 8, the comparison of errors for two types 
of nodes, by using different time steps is shown. 


Table 7: The comparison of errors of Example 3, at t = 10007 by using scattered and 


equally spaced nodes where time step is T = 10~¢. 


Equally Spaced Scattered 

N [lu — uno [e— unl lu unlnn [u—unls 

49 3.07732E — 3 6.27540E — 3 5.39392E — 3 7.59093E — 3 
64 5.49423E — 4 1.12731E — 3 1.08340F — 3 1.38534 — 3 
81 3.37897E — 4 8.97249E — 4 6.02756E — 4 1.14500E — 3 
100 5.51337E — 5 1.40970EF — 4 1.58731E — 4 2.24005E — 4 
121 3.34659E — 5 1.06018.E — 4 4.93030E — 5 1.29525E — 4 
144 5.23636E — 6 1.69188E — 5 1.28415F — 5 2.64670E — 5 
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Table 8: The comparison of errors of Example 3, by using various time steps at N = 100 


and t = 10007. 
Equally Spaced Scattered 
Time Step (F=10-)]__ Iu — unos [e—unly | leon Ta [u— uns 

© 5.51337E — 5 1.40970 —4 || 1.58731E —4 2.24005E — 4 
1/2xr 2.72150E — 5 6.72841E —5 || 6.89102EF —5 1.04333. — 4 
1/4xr 1.34159F — 6 3.27630E —5 || 3.20549F — 5 5.04380E — 5 
1/8x7 6.65022E — 6 1.61537E —5 || 1.54562B —5 2.48065E — 5 
1/16 x 7 3.30963E — 6 8.01907E — 6 || 7.58896E — 6 1.23023.E — 5 


Figure 7: The plot of the error in Example 3, at t = 10007 and N = 100 with equally 


spaced nodes. 


Figure 8: The contour plot of the error in Example 3, at t = 10007 and N = 100 with 
equally spaced nodes. 


Example 4. For the final example, we solve the following equation in a 
polygonal domain. 


0.0010 


0.0005 


-—0.0005 


-0.0010 
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v) =f, (x,y) € (-1,1)”, t> 0, 
(12) 


v(+1,y,t) = 0, v(#,+1,t) =0,t>0, 


where 7 = 100, the exact solution is v (x, y,t) = (a? — 1) (y? — 1) e-™ t and 
the function f is computed accordingly. 

In this example, for showing the flexibility of our method, we consider a 
polygonal domain with the following vertices that are arranged in clockwise 
order: 

—1)}. 


Figure 9 shows the shape of the domain and equally spaced points. As before, 
the comparison of errors for the equally spaced and the scattered nodes is 
presented in Table 9. we also plot the absolute error in Figure 10 and the 
contour plot of the error in Figure 11. 


{(—1,1), (—1,0) , (—0.286, 0), (0.5, 1), (1,1), (1, 


Table 9: The comparison of errors of Example 4, by using various time steps at N = 165 
and t = 1000r. 


Equally Spaced Scattered 
Time Step (7 = 1074) lu = Un lloo w= unlly [lu = UN lhoo ju-untlo | 
* 3.20101E — 4 1.48746E — 3 || 3.69452E —4 1.00703E — 3 
1/2xr 1.03216E — 4 4.58850E — 4 || 1.49190F — 4 3.25760B — 4 
1/4xr7 3.15830E — 5 1.28779E — 4 || 6.47536E — 5 1.11007E — 4 
1/8 x 7 9.89283E — 6 3.56560E — 5 || 3.01340E —5 4.73288E — 5 
1/16 x 7 3.35178E — 6 1.11966F —5 || 1.45773E —5 2.35022E — 5 


O.5- 


0.0 


-0.5 +] 


n n 
-1.0 -0.5 0.0 0.5 1.0 


Figure 9: Polygonal domain including 165 equally spaced points for Example 4. 


Numerical solution of Allen—Cahn equation by using ”shifted” ... 193 


Absolute Error 


mi 
mm N 


Figure 10: The plot of the error in Example 4, at t = 10007 and N = 100 with equally 
spaced nodes. 


Figure 11: The contour plot of the error in Example 4, at t = 10007 and N = 100 with 
equally spaced nodes. 


5 Conclusion 


In this paper, a fully-discrete approximation of the Allen—Cahn equation was 
considered. Based on the forward Euler/Crank—Nicolson scheme (in time) 
and the RBF collocation method (in space), we proposed a numerical scheme, 
in which the nonlinear term could be treated explicitly and the resultant 
numerical scheme was linear and easy to implement. Numerical solvability 
and stability of the method, by using second order finite difference matrices 
were discussed. Four numerical test problems were considered, which showed 
the efficiency and reliability of the proposed method. The obtained numerical 
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results showed that the equally spaced nodes are better in comparison with 
scattered nodes. 


Acknowledgements 


The authors would like to thank the anonymous referees for their valuable 
comments and suggestions, that greatly improved the final version of this 
article. 


References 


1. Allen, S.M. and Cahn, J.W. A microscopic theory for antiphase boundary 
motion and its application to antiphase domain coarsening, Acta Metall. 
27 (1979) 1085-1095. 


2. Bartels, S. Numerical methods for nonlinear partial differential equations, 
Springer Series in Computational Mathematics, vol. 47, Springer, 2015. 


3. Buhmann, M.D. Radial basis functions: Theory and implementations, 
Cambridge University Press, Cambridge, 2004. 


4. Carr, J.C., Fright, W.R. and Beatson, R.K.Surface interpolation with ra- 
dial basis functions for medical imaging, IEEE Trans. Med. Imaging 16 
(1997) 96-107. 


5. Cecil, T., Qian, J. and Osher, S$. Numerical methods for high dimensional 
Hamilton-Jacobi equations using radial basis functions, J. Comput. Phys. 
196 (2004) 327-347. 


6. Cheng, A-H.D., Golberg, M.A., Kansa, E.j. and Zammito, G. Exponential 
convergence and h-c multiquadric collocation method for partial differential 
equations, Numer. Meth. Part. Differ. Equ. 19 (2003) 571-594. 


7. Cheng, M. and Warren, J.A. An efficient algorithm for solving the phase 
field crystal model, J. Comput. Phys. 227 (2008) 6241-6248. 


8. Fasshauer, G.E. Meshfree approximation methods with MATLAB, World 
Scientific, Hackensack, 2007. 


9. Feng, X., Li, Y. and Zhang, Y.Finite element methods for the stochas- 
tics Allen—Cahn equation with gradient-type multiplicative noise, SIAM J. 
Numer. Anal. 55(1) (2017) 194-216. 


10. Franke, R. Scattered data interpolation: tests of some methods, Math. 
Comput. 38 (1982) 181-200. 


Numerical solution of Allen—Cahn equation by using ”shifted” ... 195 


11. Hon, Y.C. and Mao, X.Z. A radial basis function method for solving 
options pricing model, Financial Eng. 8 (1) (1999) 31-50. 


12. Kessler, D., Nochetto, R. and Schmidt, A. A posteriori error control 
for the Allen—Cahn problem: circumventing Gronwall’s inequality, M2AN 
Math. Model. Numer. Anal. 38 (2004) 129-142. 


13. Kim, J., Jeong, D., Yang, S. and Choi, Y. A finite difference method for a 
conservative Allen—Cahn equation on non-flat surfaces, J. Comput. Phys. 
334 (2017) 170-181. 


14. Kufner, A. and Persson, L.E. Weighted inequalities of Hardy type, World 
Scientific Publishing Co., Inc., River Edge, NJ, 2003. 


15. Lee, H. and Lee, J. A semi-analytical Fourier spectral method for the 
Allen—Cahn equation, Comput. Math. Appl. 68 (3) (2014) 174-184. 


16. Madych, W.R. and Nelson, S.A. Multivariable interpolation and condi- 
tionally positive definite functions I, Math. Comput. 54 (1990) 211-230. 


17. Mohammadi, V., Mirzaei, D. and Dehghan, M. Numerical simulation and 
error estimation of the time-dependent Allen—Cahn equation on surfaces 
with radial basis functions, J. Sci. Comput. 79 (2019) 493-516. 


18. Park, J. and Sandberg, I.W. Universal approximations using radial-basis- 
function networks, Neu. Comput. 3 (1991) 246-257. 


19. Pham, H. Handbook of engineering statistics, Springer, 2006. 


20. Shen, J., Tang, T. and Wang, L.L. Spectral methods: Algorithms, analysis 
and applications, Springer-Verlag, Berlin, 2011. 


21. Shu, C. and Ding, H. Numerical comparison of least square based 
finite difference (LSFD) and local multiquadric-differential quadrature 
(LMQDQ) methods, Comput. Math. Appl. (submitted). 


22. Shu, C., and Ding, H. and Yeo, K.S. Local radial basis function-based 
differential quadrature method and its application to solve two-dimensional 
incompressible N Navier-Stokes equations, Comput. Methods Appl. Mech. 
Eng. 192 (2003) 941-954. 


23. Sarra, S.A. A local radial basis function method for advection-diffusion- 
reaction equations on completly shaped domains, Appl. Math. Comput. 
218, 19 (2012) 9853-9865. 


24. Tolstykh, A.I., Lipavskii, M.V. and Shirobokov, D.A. High-accuracy dis- 
cretization methods for solid mechanics, Arch. Mech. 55 (2003) 531-553. 


25. Wright, G.B. and Fornberg, B. Scattered node compact finite difference- 
type formulas generated from radial basis functions, J. Comput. Phys. 212 
(2006) 99-123. 


196 C.G. Keshavarzi and F. Ghoreishi 


26. Wendland, H. Scattered data approximations, Cambridge Monographs on 
Applied Computational Mathematics, Cambridge University Press, 2005. 


27. Wu, X., Zwieten, G.J. and Zee, K.G. Stabilized second-order convex split- 
ting schemes for Cahn-Hilliard models with application to diffuse-interface 
tumor-growth models, Int. J. Numer. Methods Biomed. Eng. 30 (2014)180— 
203. 


28. Yoon, J.Spectral approximation order of radial basis function interpola- 
tion on the Sobolev space, SIAM J. Math. Anal. 33 (2001) 946-958. 


29. Yoon, J. Lp-error estimates for “shifted” surface spline interpolation on 
Sobolev Space, Math. Comput. 72 (243) (2003) 1349 —1367. 


30. Yoon, J. Improved accuracy of Lp-approximation to derivatives by radial 
basis function interpolation, Appl. Math. Comput. 161 (2005) 109-119. 


31. Zhang, J. and Du, Q. Numerical studies of discrete approximations to the 
Allen—Cahn equation in the sharp interface limit, SIAM J. Sci. Comput. 
31 (4) (2009) 3042-3063. 


